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Abstract. I present a concise review of advances realized over the past three years on planar Poisson- 
Voronoi tessellations. These encompass new analytic results, a new Monte Carlo method, and application 
to experimental data. 

PACS. PACS 02.50.-r Probability theory, stochastic processes, and statistics - PACS 45.70.Qj Pattern 
formation - PACS 87.18.-h Multicellular phenomena 
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In this talk I will concisely review a coherent collection of 
new results on planar Voronoi cells obtained over the last 
three years PQ[H[31I3]. In Fig. [1] you see a Voronoi tessella- 
tion: the set of dots is given and I will refer to them as 
'centers' or 'point particles.' Perpendicular bisectors have 
been constructed to the line segments connecting nearby 
particles. The bisectors meet in trivalent vertices (unless 
there is accidental degeneracy) and partition the plane 

into convex polygons called Voronoi cells. tures, whether occurring naturally or synthesized; exam- 

In physics Voronoi cells have, roughly, two broad classes pies are biological tissue or soap froths, (ii) They may 
of applications, (i) They may directly model cellular struc- serve as a tool of analysis. For example, in a dense con- 



Fig. 1. Voronoi tessellation of a set of point particles. 
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figuration of real physical particles they define lattice de- 
fects; and if you decide to study your favorite theory on 
a lattice of randomly located spatial sites, the Voronoi 
construction is the natural way to define nearest-neighbor 
relations between these sites. Longer lists of applications Q 5 10 15/2 

may be found in the literature [SKHIIS]- 



When the point particles are distributed independently 
and uniformly - which I will henceforth stipulate the 
resulting tessellation is called a foisson- Voronoi tessella- 
tion. Real systems will of course most often exhibit devi- 
ations from this simple mathematical model: soap froths 
evolve with time and their cells are not constructed around 
a center; real particles have a minimum distance of ap- 
proach; etc. However, the Poisson- Voronoi tessellation, be- 
ing the simplest one to study, is to cellular systems what 
the Ising model is to magnetism. 



Fig. 2. Histogram of the sidedness probability p„. 

expression for p„ is on the same footing as a partition 
function of n interacting particles in two dimensions. No 
exact evaluation is known, but Monte Carlo simulations 
have led to the histogram of Fig. [21 It shows that p„ peaks 
at rt = 6 and decreases very rapidly for larger n. It is as 
small as 10~^ for n ~ 15 and 10^^" for n sa 18. Drouffe 
and Itzykson (DI) [8, 9 devised a Monte Carlo algorithm 
that yielded Pn with a single-digit accuracy up to n 25. 
The asymptotic decay of p„ for large n has been a subject 
of speculations. Many authors fitted p„ to a decaying ex- 
ponential or stretched exponential, often multiplied by a 
power. DI proposed p„ ~ n^"" with a w 2. 



Over the last 50 years, starting with the work of Meij- 
ering [7|, many analytical properties of the Voronoi cell 
have been determined. These include the statistics of the 
perimeter segments, of the angles at the vertices, the cell 
area, and so on (an overview is given in Ref. [5]). How- 
ever, the Voronoi cell's most prominent property, viz. its A new line of research began with the mathematical chal- 



2 Mathematics 



probability p„ of being n-sided, is extremely difficult to lenge of obtaining the analytic asymptotic expression of 

access analytically. The basic expression for p„ is a 2n- pn for n ~* oo. The calculation, although involving only 

dimensional integral on the positions Rm {m = 1, . . . , n) classical tools of analysis, is of considerable complexity [J 
of n point particles neighboring a central one placed in I present briefly its main ideas. 

the origin. A coupling between the Rm arises from the A key role is played by the angles and rji defined 

condition that the resulting cell be n-sided. Hence the in Fig.[3l They define the cell up to a radial scale factor. 
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After writing p„ in terms of the variables of integration 
^rn. and r]i, one can express it in the factorized form 

V -C - ^^^""^^""^ (2) 

Pn — '-^nPn i Pn — (^2n)! ' 



where pn^ results from integration with respect to the P^^^ 
and P2^'^ . Hence the problem of finding p„ has been re- 
placed with that of finding C„ . This problem involves the 
Fig. 3. The angles ^„ and r^e (with m,f = 1, . . . ,n) that define 0{n~^/^) and 0{n-^) corrections to the infinite-n behav- 
the n-sided Voronoi cell up to a radial scale factor. ior. It appears that for C„ one can set up a perturbation 

expansion which shows that 

The obvious sum rules X]m=i ^™ = ^tt and Vt = 27r 

C„ = C[l + O(n-i)], n^oo, (3) 

represent infinitely weak constraints in the limit n — > oo. 

It appears that in that limit, to zeroth order in n^^, the where C can be expressed elegantly as 

angles and ry£ become a set of 2n independent variables, °° / 1 4 \ 

(0) ro) ^=11(1-72+7) = 0.344347... (4) 

the distributions Pf (n^„i/27r) and P2 {nrie/2TT) being q=i ^ ^ 

given by This infinite product is in fact the 'elastic' configurational 

partition function of the cell perimeter. The factor of index 
^1 = 4a;e \y) — ^ ^ ^ x,y > 0. (1) q in (jl]) stems from deviations from circularity having a 

wavelength 2'KRc/q along the perimeter. 



This independence is nontrivial: the S,m and rj^ overlap, 
so how can they be uncorrelated? This property is best 
initially introduced as a hypothesis, to be confirmed self- 
consistently later by the solution procedure. 

The independence of the angles can be shown to imply 



There is an interesting corollary. It says that there is 
a continuum limit in which the cell perimeter, expressed 
as a function P(0) of the polar angle satisfies 

, ,2 = L{(t>), (5) 



that as n — > 00 the shape of the n-sided cell approaches where L[(j)) is Gaussian noise which is white to order zero 

a circle with probability 1. From an analysis of the radial but has colored order corrections. The solution of (O, 

integral in the expression for p„ one concludes that (for under appropriate conditions pertaining to the average of 

point particle density A) this circle is of the radius Rc{n) = ^(0) and to its periodicity [5], produces the full functional 

-\/ n/AirX. probability distribution V[R{(I))] of an arbitrary perimeter 
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R{4>). Hence the asymptotic determination of p„ leads to the form 

a complete understanding of the behavior of the large n- Cn = {WnOn), (6) 

sided cell. 



Finally, the mathematical methods employed here are Cqj (Q^ 

respect to the and P2 ■ Furthermore Wn is a weight 



valid for all n = 3,4, . . .. Here (. . .) is an average with 
applicable to other problems that have arisen in mathe- 

that can be expressed [1], through a series of equations, 

matics and in physics. I mention two of them. 

in terms of the sets of angles {^m} and {77^}, and On is a 

The Crofton problem. Let a plane be traversed by in- ^ 

projector: &„ = 1 if a certain condition on the angles is 

tersecting straight lines, distributed randomly and uni- 

satisfied, and &„ = if not. This condition requires the n 

formly. This is another way of partitioning it into convex 

points Rm all to be located such that they contribute a 

cells. The 'Crofton problem' is the question of determin- 

nonzero segment to the perimeter. All previous 'fixed n' 

ing the probability p^'^ for the cell containing the origin to 

Monte Carlo methods described in the literature run into 

be n-sided. This was done in Ref . [IQ4 1 again in the limit 

such an acceptance criterion. In all proposed cases, the 

of large n. 

projector rejects an ever larger fraction of configurations 

The problem of extremal points. Let n points be dis- 

as n gets larger: this is the phenomenon of attrition, well- 

tributed randomly and uniformly in the unit disk. Then 

known (and deplored) in many Monte Carlo studies. 

can one determine the probability p'^^ that the convex 

The particular split ^ that I made into a zeroth order 

envelope of this set is an n-sided polygon? Work on this 

problem and a remainder now appears to be the right one 

problem is in progress. 

[Ij. If the £,m and rjg are chosen from the zeroth order 
distributions H]), then the problem of attrition in the large- 

3 Monte Carlo method , , , , r ■ r 

n limit is eliminated: as n grows, the fraction of accepted 

Attempts to estimate p„ by Monte Carlo simulation have configurations tends to unity! In practice, less than 1% of 

been numerous (see Ref. fSl). Straightforward methods star- the configurations is rejected for n > 20, and less than 

ting from a random set of centers suffer from the fact that 0.01% for n > 40. 

many-sided cells are very rare. More sophisticated meth- Two things then become possible, 

ods like those of Ref. [8J aim at directly generating n-sided (i) To accurately determine the sidedness probabili- 

cells for an n fixed in advance. ties pn for arbitrarily large n. Tables with values of p„ 

The derivation of the asymptotic result ©-([l]) for C„ to at least four decimal places are given in Ref. [1] in the 

has as its starting point a non asymptotic expression of range 3 < n < 1600. For high n these probabilities be- 
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io[ 



come unphysically small: one has pioo — 5.269 x lO^^^'* 
and piooo = 6.36 x io~384i (^j^/)^ j-^^ j gj^ow below 
which rewards can be gained from studying them. 



(ii) To Monte Carlo generate typical n-sided cells for 
arbitrary a priori given n, together with their 'natural' 
environment of other cells. This is done as follows. An- 
gles and rji are randomly picked from the zeroth order 
distribution (but subject to the sum rules). If they pass 
the acceptance criterion, the cell perimeter is constructed. 
By radial scaling the cell radius is given its most proba- 
ble value Re- The positions of the central point particle 
and its first neighbors are now fixed. Next, all other point 
particles, in order that they do not interfere with the cen- 
tral Voronoi cell, are excluded from the 'fundamental do- 
main' associated with that cell, i.e. from the union of the 
n disks centered at the cell vertices and passing through 
the origin; however, those other point particles (which are 
second and higher order neighbors to the central one), 
occupy the remaining region of the plane uniformly and 
randomly. Hence, after Monte Carlo generation of their 
positions, the construction of the full Voronoi tessellation 
may be completed by any of the existing algorithms. 



Fig-Sshows a 96-sided Voronoi cell constructed in this 
way. Snapshots of cells with n = 24, 48, 96, and 1536 may 
be found in Ref. [4]. 




-10 10 

Fig. 4. Snapshot of a 96-sided Voronoi cell. 

4 Aboav's law 

Aboav's law, formulated in 1970 [11], holds that the neigh- 
bor of an n-sided cell has itself an average number m„ of 
sides given by 

m„ = a + bn~^, (7) 

where a and b are positive constants. Eq. ([7]) used to be, 
and is still often, formulated as nm„ = an+b, and is there- 
fore called the 'linear law'. Since m„ decreases with n, it 
says that "many-sided cells tend to have few-sided neigh- 
bors, and conversely. " Two-cell correlations have been mea- 
sured in a large variety of experimental systems and found 
in agreement with Aboav's law for parameters in the range 
4-5 ;$ a ^ 5.3 (in this context 5 has been called a 'magic 
number') and 5.7 < 6 < 8.5. 

The first one to present heuristic arguments for the 
validity of ([7]) was Weaire [T2], whence the alternative 
name 'Aboav- Weaire law.' Other 'proofs' of this law (see 
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Ref. [3] for some references) appeal to mean field approxi- 
mations or proceed by maximizing a hypothesized entropy 
functional. While some workers consider Aboav's law as ^ 

r 

no more than a linear approximation to some unknown 
curve, others have attributed to it a more fundamental 
status, as is clear from numerous statements in the litera- 
ture: "In this paper we present an alternative derivation of 
the Ahoav-Weaire law" |13| . or "In all known natural and 
man-made structures, it is found empirically that nnin is 
linearly related to n" (14j, and many others. 



Brakke (1985) simulation 
Aboav's law: m = a + b/n 



This work: m^ = 4 + 3 (n/n) 

1/n 



1/2 



0.3 



Fig. 5. Monte Carlo data, Aboav's law, and this work's asymp- 
totic expression ([8| for m„. 



Now it has been known for over twenty years that 
Aboav's law does not hold exactly for the Poisson- Voronoi 
tessellation. This is borne out by Fig.O numerical simu- 
lation shows unambiguously that the m„ curve has a very 
small but distinct downward curvature. It is so small that, 
if also present in nature, most experimental data would Fig. 6. Zoom onto the box in Fig.|l showing some first- 




not detect it. 



In face of this, the defenders of Aboav's law as a ba- 
sic truth hold that the Poisson- Voronoi tessellation, be- 
cause of being constructed around random centers, does 
not correspond to any true physical structure. In real cel- 
lular structures, so they argue, any initial randomness has 
always undergone a relaxation process (for example to re- 
lieve internal stress); and this then would have induced 
the validity of Aboav's law. 



neighbor cells to the central cell. 

4.1 Aboav's law and the Poisson-Voronoi tessellation 

Before returning to real physical systems, let us first see 
what the present theory, further developed in Ref. [3] , has 
to say about the validity of Aboav's law for the Poisson- 
Voronoi tessellation. The picture of Fig. [6] is useful for the 
argument. It shows a phenomenon that becomes ever more 
pronounced as n gets larger (whence the need to consider 
very large n; see 0] for a picture with n = 1536): the first 
neighbors of the central cell become more and more elon- 
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gated while their width goes down to order 1 / y^. Since 
the second neighbors of the central cell have dimensions of 
order unity (for unit particle density), the geometry dic- 
tates that most first neighbors must become four-sided. 
The fraction of them that is five-sided is only of order 
1 / this is because five-sidedness occurs only when a 
first-order neighbor is adjacent to two second-order neigh- 
bors. In Fig. [6] all five-sided cells have been marked. 

It follows that m„ is equal to 4 plus 0{n^^^'^) correc- 
tions. The latter require a nontrivial calculation [3] and 
one finds 

m„ — A + ^^/tt/u + . . . , n — > oo. (8) 

In Fig.Othe Monte Carlo data for rn„ (due to Brakke [T5] ) 
are compared to the first two terms of the asymptotic ex- 
pansion ((8]), as well as to the best linear two-parameter 
fit of type ([7]) . Clearly the asymptotic theory that I devel- 
oped explains for the first time why m„ must be curved; 
moreover, it predicts correctly the order of magnitude of 
this curvature [3]. 

4.2 Aboav's law and experimental systems 

Under suitable conditions, polystyrene latex spheres of di- 
ameter « l/iTO, when trapped at a water/air interface, 
will on a time scale of hours undergo a process of dif- 
fusion limited colloidal aggregation (DLCA). The exper- 
iment was performed by Earnshaw and coworkers [161 
[TTifTS] in the 1990's and was simulated very recently by 
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Earnshaw and Robinson 




1/n 

Fig. 7. Experimental and simulational m„ data compared to 
the asymptotic Poisson- Voronoi curve ((Sj. 

Fernandez- Toledano et al. [19^. Snapshots of the system 
were taken at regular time intervals and the center of grav- 
ity of each cluster of particles was determined. Then the 
Voronoi cells belonging to this set of centers were con- 
structed and m„ was obtained. Fig. [7] shows the experi- 
mental DLCA data (dots) of Ref. [TO] taken after 60 sec- 
onds, as well as the corresponding simulation data (di- 
amonds) of Ref. [in]. Also shown is the same theoretical 
curve (solid line) as in Fig. [5] In view of the experimen- 
tal and numerical error bars, this figure does not prove 
that the DLCA system is of the Poisson- Voronoi type. I 
believe however that, as an explanation for these data, 
the present zero-parameter curvature-predicting theory is 
preferable to the linear two-parameter fit ([7]). 

A final point is worth discussion. In the experiment, 
the clusters whose centers of gravity serve as the 'point 
particles' of the Voronoi construction, have some finite 
diameter d. Is this a problem? The experiment is done. 
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of the present asymptotic theory can be expected to be 
applicable only for n < n* , where n* is a sidedness such 
that w d. The point particle/cluster density A drops 
out of this relation and we find n* 100. Taking into 
account the possibly fractal structure of the clusters may 
lower this estimate of n* , but the n values of Fig. [7] are less 
than n* by what seems a safe margin. Crossover between 
the ideal Poisson- Voronoi behavior and a finite-diameter 
theory was briefly discussed in Ref. [3]- 



What we learn from this is that in other experimental 
situations, too, as the range and the precision of the data 
increase, a curved rn„ law may well be expected and must 
indeed be looked for. 
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